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Abstract. We analyze a one dimensional quantum system consisting of a test particle inter- 
acting with two harmonic oscillators placed at the positions a±, 0,2, with a\ > 0, \a^\ > ax, in 
the two possible situations: 02 > and a-i < 0. At time zero the harmonic oscillators are in 
their ground state and the test particle is in a superposition state of two wave packets centered 
in the origin with opposite mean momentum. Under suitable assumptions on the physical pa- 
rameters of the model, we consider the time evolution of the wave function and we compute 
the probability V~ lTl2 {t) (resp. V^ in2 (t)) that both oscillators are in the excited states labelled 
by ni, tt-2 > at time t > \a2\vQ 1 when 02 < (resp. 02 > 0). 

We prove that V^ in2 (t) is negligible with respect to V£ Ln2 (t), up to second order in time 
dependent perturbation theory. 

The system we consider is a simplified, one dimensional version of the original model of a cloud 
chamber introduced by Mott in [M], where the result was argued using euristic arguments 
in the framework of the time independent perturbation theory for the stationary Schrodinger 
equation. 

The method of the proof is entirely elementary and it is essentially based on a stationary phase 
argument. We also remark that all the computations refer to the Schrodinger equation for the 
three-particle system, with no reference to the wave packet collapse postulate. 



1. Introduction, notation and result 



In his paper of 1929 Mott QMJ) analyzes the dynamics of formation of tracks left an a-particle 
emitted by a radioactive source inside the supersaturated vapour in a cloud chamber. He 
notices the difficulty to understand intuitively how a spherical wave function, describing the 
particle isotropically emitted by the source, might manifest itself as a straight track in the cloud 
chamber. 

Without referring to any wave packet collapse, he proposes an explanation based on the analysis 
of the whole quantum system made up of the a-particle and of the atoms of the vapour. Using 
a simplified model with only two atoms and making use of time independent perturbation 
arguments, he concludes that each ionization process focuses the probability of presence of the 
a-particle on narrower and narrower cones, around the straight line connecting the source to 
the ionized atoms. 
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In this way Mott suggests a quantum dynamical mechanism responsible of the transition be- 
tween an initial superposition of outgoing waves heading isotropically in all directions toward 
an incoherent (classical) sum of those same waves. 

We mention that the same problem is also discussed in [H] and later in |Bej . where the above 
approach is compared with the explanation based on the wave packet collapse. We refer to |Brj . 
|HA] . [CL] . [BPT] for some further elaborations on the subject and to |LR] for a description of 
the original experimental apparatus. 

The aim of our work is to provide a detailed time dependent analysis of a one dimensional 
version of the system investigated by Mott. The system we consider consists of a test particle 
and two harmonic oscillators. In our model a superposition of two wave packets centered in 
the origin with opposite momentum plays the role of the spherical wave of the a-particle and 
the oscillators replace the atoms to be ionized. Under suitable assumptions on the physical 
parameters of the model, we perform a detailed time analysis of the evolution of the system 
wave function using time dependent perturbation theory and we give a quantitative estimate 
of the joint excitation probability of the oscillators. Roughly speaking, our main result is that 
such probability is essentially zero if the oscillators are placed on opposite sides of the origin, 
while it has a finite, non-zero value in the other case. Following the line of reasoning of Mott, 
the result can be interpreted saying that before the interaction the test particle is delocalized 
while after the interaction it is either on the left (if there is an excited oscillator on the left) 
or on the right (if there is an excited oscillator on the right). In any case one can say that the 
test particle propagates along an almost classical trajectory, without making any reference to 
the wave packet collapse postulate. 

In |CCF] the authors consider a similar problem in three dimensions where a particle interacts 
via zero range forces with localized two level quantum systems. A non perturbative analysis of 
the model is carried out but results are valid only in the scattering regime. 
Let us introduce the model. We consider a three-particle non relativistic quantum system in 
dimension one, made of one test particle with mass M interacting with two harmonic oscillators 
with the identical mass m. We denote by R the position coordinate of the test particle and by 
r\,T2 the position coordinates of the two oscillators. The Hamiltonian of the system in L 2 (IR 3 ) 
is 

H = H + XH 1 (1.1) 

ft 2 h 2 1 ft 2 l 

H = A R A n + -mw 2 (r x - ax) 2 A r „ + -mu 2 (r 2 - a 2 ) 2 (1.2) 

2M 2m 1 2 v ; 2m 2 2 v ; v ; 



H x = V(S- l (R - n)) + V(S~ l (R - r 2 )) (1.3) 

where A > 0, u > 0, ai > 0, a 2 G JR., with a\ < \a 2 \, 6 > and V is a smooth interaction 
potential. The assumptions on V will guarantee that the Hamiltonian H is self-adjoint with 
the same domain of Hq and then the evolution problem corresponding to the Hamiltonian H 
is well posed. For the test particle we choose an initial state ifj in the form of a superposition 
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state 



ifj(R) = i(j + (R)+^-(R) (1.4) 
^(R) = ^=e-^e ±l ^ R , P = Mv (1.5) 



jo 

' Pn<?\ 2 



where o > 0, Af = [2-^(1 + e~{~*~) )] 1//2 is the normalization factor and Pq,vq denote the 
absolute value of the initial mean momentum and velocity of the test particle. 
For the harmonic oscillator centered in a,j, j = 1,2, the initial state (f>Q J is the corresponding 
ground state. Moreover we define 



^ nj {r\r 3 -a 3 )) (1.6) 



7=a/— (1-7) 
V mijj 

where m is the normalized Hermite function of order m G N. We notice that the parameter 7 
has the dimension of a length and gives a measure of the spatial localization of the oscillators. 
Let us denote by n, r 2 , t) the wave function of the system; ty(t) = \l/ (•,-,-, i) is the solution 
of the Cauchy problem 

^^(*) = H*{t) (1.8) 



v&(O)=^0S 1 ^ 2 (1.9) 

We are interested in the probability that both harmonic oscillators are in an excited state at 
a given time t > 0. The solution of the three-body problem fll.8p . fll.9l) is not known in closed 
form; we shall limit ourselves to a perturbative computation. It is worth mentioning that, in 
order to get a non trivial result, we are forced to compute the second order approximation of 
the solution of the Cauchy problem fll.8p . (11.91) . which we denote by ^(-R, fi, r 2, t). Therefore 
the object of our analysis is the quantity 



dR 



dndr 2 < (n)#* (r 2 )* 2 (i?, r 1; r 2 , t) 



;i.io) 



for ri\ 7^ and n 2 7^ 0, where ± refers to the cases a 2 > and a 2 < respectively. Formula 
fll.lOP represents the probability that both oscillators are in an excited state at time t, up to 
second order in perturbation theory. 

The explicit computation of fll.lOp will be performed exploiting some further assumptions on 
the physical parameters of the model. More precisely the complete set of assumptions required 
for our analysis is the following 
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(Ao) 



A ^«l (1.11) 



(Ax] 

quai 
all 0(e) where 



The quantities Sm=§, 8E = 6R= fa 6L = ^ 6r s = jfa, for j = 1, 2, are 

A < e < 1 (1.12) 



The interaction potential V : R — > R a continuous, positive and compactly supported function. 

Let us briefly comment on the above assumptions. In (A Q ) we ensure that the dimensionless 
coupling constant Ao is small. In (A\) we assume that the mass and the kinetic energy of the 
test particle are much larger than the mass and the spacing of the energy levels of the oscillators; 
moreover the initial wave packets of the test particle are assumed to be well localized and the 
interaction is required to be short range; finally the characteristic time of the oscillators u~ l is 
assumed to be much smaller than the flight times r±, r 2 of the test particle, which are defined 
by 

°1 l a 2| /-, 1Q \ 

n = — , r 2 = (1.13) 

Condition ( 11. 12ft guarantees that the first and second order corrections in perturbation theory 
remain small compared with the unperturbed wave function, in fact of order Ao£ _1 and \q£~ 2 
respectively. 

In order to understand the meaning of (Ai), let us consider the parameters M, v , a\ a 2 all of 
order one. Then we obtain m = 0(e), u = 0(e~ l ), hu = 0(e), o = 0(e), 5 = 0(e). 
We observe that the length 7 introduced in (II. 7p can be written as 



7 = ai(W^ (1.14) 

and this means that 7 is of the same order of 6 and a. In particular this guarantees that the 
transit time of the test particle on the region where the oscillators are localized is of the same 
order of the characteristic time of the oscillators. To simplify the notation, from now on we 
shall fix 

5 = 7 (1.15) 

We also introduce here a (large) parameter which is useful to express the various estimates in 
the proof 
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A^^ = 0(e- 1 ), j = 1,2 (1.16) 
Our main result is the following. 

Theorem 1. Lei us assume (A ), (Ax), (A 2 ) and fix t > r 2; nx ^ 0, n 2 ^ 0. T/ien /or any 
fcGN, wift k > 2, we have 

l /V 4 



P»it*(*)< Xf^(f) C -n 2 (*) (1-17) 



^ 2 W = 16vr 4 v^(^)V 2 



11 %i)(^%>o)(gj 

J'=l,2 



2 

+ 5 nin2 (t) (1.18) 



= (1.19) 

i Ao v4 



|5 ni „ 2 (t)|< — 1-^1 D nin2 (t) (1.20) 

where the symbol ~ denotes Fourier transform and Cnimif), D niTl2 (t) are functions of the 
physical parameters of the model which will be explicitely given during the proof (see 
qjjty below). 

We remark that the estimates (11.171) . (11. 18ft . ( 11.201) are not optimal; in particular Cnin 2 {t) and 
D ni n 2 (t) diverge for t — > oo. From (I4.16p . (14. 19j) it will be clear that Cnin 2 (t), D nin2 (t) are of 
order one, and then the estimates are meaningful only for t larger but of the same order of r 2 . 
Let us briefly outline the strategy of the proof and give a heuristic argument which, at least at 
a qualitative level, justifies the result stated in theorem 1. We find convenient to represent the 
solution of (II. 8p . (II. 9p in the form 

*(R,n,r 2i t) = f^n 2 (R,t)4>n\(ri)4>nl(r2) (1-21) 

ni,n 2 

where f ni n 2 (-,t) — fn in2 (t) belongs to L 2 (M) for any ni,n 2 G N and t > 0, and it is explicitely 
given by 



f nin2 (R,t) = J drxdr^irxW^i^MR^r^t) (1.22) 

We notice that the coefficients of the expansion f nin2 (R,t) have a precise physical meaning; in 
fact the quantity 

r dR\f nin2 (R,t)\ 2 (1.23) 
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represents the probability to find at time t the test particle in Q C R when the two oscillators 
are in the states labelled by nx,n 2 respectively 

The equation for the coefficients f nin2 (t) is obtained from Duhamel's formula 

Wt) = e-^ ^ - *~ / dse-T^-^H^h) (1.24) 
n Jo 

multiplying by (p^^nl an d then integrating with respect to the coordinates of the oscillators. 
The result is 

fn im (t) = fW n2 (t) - JjsT nin2 {t - s) ^ V:i n f nn2 {s) + £ V:* j2 fn lh (s)J (1.25) 
where in the above formula we have introduced the notation 

/SJ*) ^^o^oe-^e-^V (1-26) 

h 2 

K = — A R (1.27) 

T nin Jt) = e -frK (L2g) 

n 



V%Jx) = Jdy^{y)^{y)V{r L {x-y)), m,neN, z = l,2 (1.29) 

We want to give an estimate of the solution of fl 1 . 2 5 [) up to second order in perturbation theory. 
Iterating twice equation HI . 25[) we obtain 



fmnM = fi% 2 (t) + f$ n2 (t) + fi% 2 (t) + E nin2 (t) (1.30) 



where 



/£!»(*) = -J/ s - S ) l^r v:^i( S ) + y: v%jw 2 (s)j (1.31) 

& (*) = - [dsT niri2 {t - s) K)JnU s ) + E V ZJZM j (1-32) 
and £ nin2 {t) is the error term which we shall neglect in the sequel. Obviously we have 



ni,n2 



Exploiting the explicit expression dL26D of /&(*), we can write /nin 2 (0 i n the form 
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fnin2 (^) 



-5 n20 I ds r ni0 (t - s) Co e~ 2 ^ e-^^V - « 



mo 



ds T 0n2 (t-s)V^ 



(1) 



(1.34) 



From formula (11.341) it is clear that fn±n 2 {t) = if n\ ^ and n 2 ^ 0. As expected, this means 
that the probability that both oscillators are in an excited state is zero up to first order in 
perturbation theory. As a consequence, from (11.101) we get 

Km(t) = JdR\fi% 2 (R,t)\ 2 , m^0,7i2^0 (1.35) 

Following the original strategy of Mott, a crucial point of the analysis is the explicit evaluation 
of fnjf,(t) and fj^ (t). We notice that V^ (x) and (erfc^tp^fe) are essentially different 
from zero only for x ~ d\ and x ~ ±VqS respectively. This means that the only non zero 
contribution to the time integral defining /„*o(t) comes from and such contribution is 

essentially concentrated around s = ^ = t x . Hence we can argue that fj^ Q (t) is approximately 
given by a wave packet starting at time T\ from the position a\ of the first oscillator, with a 
velocity close to Vq. In particular it is essentially different from zero only in a neighborhood of 
a>i + vo(t — Tx), for t > T\. 

Analogously, f^ (t) is approximately given by a wave packet starting at time r 2 from the 
position a 2 of the second oscillator, with a velocity close to v if a 2 > 0, and to — v if a 2 < 0. 
Then fj^ (t) is essentially different from zero only in a neighborhood of a 2 + v (t — r 2 ), for 
t > r 2 , a 2 > 0, and in a neighborhood of a 2 — Vo(t — r 2 ), for t > r 2 , a 2 < 0. 
Let us now consider the second order term /?lin 2 (t); exploiting expression (ll.34p . we have 



Jo „■ Jo 



Jl J2 



dsT nin2 (t - 8)Y% f®(8) - fdsT nin2 {t - s)V:* fi%(s) (1.36) 

Jo 

Since we are interested in the probability that both oscillators are excited, only the last two 
terms of (11.361) are relevant. 

We notice that the supports of V^ and f^Joi 8 ) are essentially disjoint for any s > and this 
implies that the third term in the r.h.s. of fl 1 . 3 6 [) gives a negligible contribution. 
For the same reason, the fourth term in the r.h.s. of (11.361) is also approximately zero if a 2 < 0. 
This explains why we expect that an estimate like (11.171) holds. 
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On the other hand, in the case 02 > the product V^ofn io( s ) i s different from zero for s ~ r 2 . 
In such case the fourth term in the r.h.s. of H 1 . 3 6 1) gives a non zero contribution and this 
explains why we can expect that a formula like fll . 18j) holds. 
We collect here some further notation which will be used in the paper. 

- (x) denotes (1 + 2 2 ) 1 / 2 ; 

- d\J is the derivative of order k with respect to x\ of a smooth function f(xx, . . . , x n ), for 
nGN and I — 1, ... ,n; 

- 11/11**.* = Ei =1 Z k m= ofdx(x) s \Kf)(x)\, keN,s>0; 

- Wfhl = \\f\\ w ^, 

- c is a generic positive numerical constant. 

The paper is organized as follows. In section 2 we study the first order approximation step. 
In section 3 we analyze the second order approximation, distinguishing the two case a 2 > 
and (22 < 0. In section 4 we compute the joint excitation probability of the two oscillators 
concluding the proof of theorem 1. Finally in the appendix we give a proof of a technical 
lemma. 

2. First order approximation 

In this section we fix t > Tj, j = 1, 2, and we give an estimate of the first order terms fn\{t). 
We only give the details for the case ai > since the opposite case can be treated similarly. 
We rewrite /^(t) as follows 

C(t) = fiT(t) + fiX(t) (2.1) 

= -r % o(t) fds e m ^ s e***> V ^ o e-i^o^i (2 . 2) 
Moreover let us define for j — 1,2 and s, t > 

I j (8)=ei' K °V% e-i' K ° (2.3) 

hf(t) = I ' dsjwijis)^ (2.4) 
Jo 

As a first step the operator (12.31) will be written in a more convenient form. 
Lemma 2.1. For any f G L 2 (M) and s > the following identity holds 

(Ij(8)f)(R) = U gA0f(R+ (M^hsOe^ e 1 ^ e~ iA ^ (2.5) 



where 
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(2.6) 



Proof. Exploiting the explicit expression of the free propagator we have 



U l dx V^{x) Idy f(y) e^v^v 



2nhs J 
■ U e~^ R2 [dyf(y)e^y 2 jdxV^x) e<^-^> 



2nhs 

= ~^- s e~^ R2 Jdy f(y) \C^M{hs)-\y - R) 

= ^ jdi f(R+ {MiT^hs i) e^ +l ^ e Vg(7 _1 (2-7) 

where in the last line we have introduced the new integration variable £ = Mj(hs)^ 1 (y — R). 
Using the convolution property of Fourier transform we have 

v n-o( x ) = - I d y (pn 3 {i~\y - s)) ^0(7^(2/ ~ a i)) v r (7 _1 - y)) 



7 

dz V(z) (0 nj o )(7~ 1 (a ; - - jz)) 

= 7 Jdk Vfrk) (K^)ink) e<*-^ h (2.8) 
Hence 

V% (k) = V^l V(lk) (4^o)(lk) e~^ k (2.9) 
Using (12. 9p in (12.71) and introducing the large parameter Aj = — we conclude the proof. 

□ 

Using the above lemma we can rewrite also the integral in (12. 4ft . 
Lemma 2.2. 



hf(t) = fds ld£ Ff(;Z,a) e iA ^,s) ( 2 .10) 
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where 



Ff(i2,C,s) = ^(0eW s 1>t{R,S,8) (2.11) 

^i2,e,«) = 4=e"^- ±,+fl (2-12) 
\/cr 



A! = — ^, A^^oi-e (2.13) 

M7 7 



#fc*) = (±f (2-14) 

and qj has been defined in 

Proof. The proof is trivial if we notice that 

ip ± (R+ (M^)- l hs^) = 4>±{R,£,s) e ±iA ^i (2.15) 

and use lemma 2.1. 



□ 



The next step is to estimate (I2.10p . i.e. an integral containing the rapidly oscillating phase 
Aj9^(^, s). The standard stationary (or non-stationary) phase methods can be used to obtain 
the estimate. 

It is worth mentioning that the integral in (I2.10p contains also other phase factors depending 
on (£, s) which, however, are slowly varying under our assumptions on the physical parameters 
of the model. 

The asymptotic analysis for Aj — > oo is simplified by the fact that 9j(£,s) is a quadratic 
function. The only critical points of the phase are (±t/j, ±Tj) and, moreover, the hessian matrix 
is non degenerate, with eigenvalues ±7"^ . This means that the behaviour of (I2.10p for Aj — > oo 
in the case with 6j is radically different from the case with due to the fact that in the first 
case the critical point never belongs to the domain of integration while in the second case this 
happens for t > Tj. 

For the analysis of this second case it will be useful the following elementary lemma. For the 
convenience of the reader a proof of the lemma will be given in the appendix. 

Lemma 2.3. Let us consider for any A > 

.7(A) = Jdx J'dy f{x, y) e iA *« (2.16) 
where fx, v are positive parameters, f is a complex-valued, sufficiently smooth function. Then 
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J(A) = ±£ X (A) (2.17) 

= I 2 vr/(0,0) + ^/C 2 (A) (2.18) 

= i27T/(0,0) + -^2«4rf y /(0,0) + ^/C 3 (A) (2.19) 
where /Q(A), / = 1,2,3, are explicitely given (see the appendix) and satisfy the estimates 

|£i(A)| < a (\\K,0)\\v + Jdx\\d x d y f(x,)\\ L )j (2.20) 
|/C 2 (A)| < c 2 (||^/(-,0)|| L1 + ||4^/(-,0)|| xl + jdx \\dldlf(x,-)\\ L ^ (2.21) 



and the constants ci, c 2 , C3 depend only on fi, v. 



Exploiting lemma 2.3 we obtain the following asymptotic behaviour of ( 12.101) for t > Tj when 
the phase is Of. 



Proposition 2.4. For any t > Tj we have 



9ttt- 1 
KV) = -j± e~ iAj9j Fj~(- , qj,^) + ^nf(; t, A,) (2.23) 



where 

^ \dlF+(R, C, Tj)\ + fdi \d^d s F+(R, C r,-)| 

(2.24) 



|7e+(i?,t,A,)| <c,- 



h /de( ids\djd 2 s F;(R,z, s )\ 2 



and Cj depends on t and r, 



j ■ 



Proof. Let us introduce the change of coordinates x — £ — qj, y = s — Tj in (I2.10p and the 
shorthand notation F(x, y) = e~ tAjq: > Ff(R, x + qj,y + Tj). Then 
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t-Ti 



h+(t)= dx dyF(x,y)e T i y (2.25) 



The integral in (12.251) has the same form as the integral (I2.16P analysed in lemma 2.3, if we 
identify u, //, /, A with Tj,t — Tj, F, — respectively. Then, exploiting formula (12.181) . we obtain 
the r.h.s. of fl2T23l with 

UUR, t, Aj) = -rf ( dx n*,0)-F(0,0)-<lF(*,0)* / + 



X 2 \ t Tj Tj 



X 



rfl j, ~ W) - <W»#» (2 . 26) 



y 

Using (I2.2ip we immediately get the estimate (12.241) and this concludes the proof. 



□ 



In the next proposition we shall analyze the asymptotic behaviour of (12.101) when the phase is 
6j. Taking into account the error term in (I2.23p . it is sufficient to show that hj(t) = 0(Aj 2 ); 
on the other hand we remark that, following the same line, it is easy to extend the result to 
hj(t) = 0(Aj k ), for any integer k. 

Proposition 2.5 . For any t > we have 



AJ 3 



where 



(2.27) 



\Kj(R,t,A3)\ < jds jdt\d\Fr(R^, s )\ 



(2.28) 



Proof. If we notice that 



-iAj ' * 



d:< ■ 



2 



(2.29) 
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and integrate by parts two times in the r.h.s. of (12.101) we easily obtain the r.h.s. of ( 12. 27ft with 



Kj(R,t,Aj) = -rf [ ds 1 / (d\F-{R^ lS )) (2.30) 

JO [S + Tj) J 

Then by a trivial estimate we conclude the proof. 

□ 



Collecting together the results of propositions 2.4 and 2.5 we finally obtain an asymptotic 
expression for t > Tj of the first order terms when Aj — > oo 

A {1) 1 

/$(*) = i~ e-^t + (2.31) 

= _27Tz ^ e -K+iM-*A^+^ g | ( } (2 32) 



^ + = ^+(-,g J ,r i ) (2.33) 
^fO.t.A,-) = -r ri . (t) {Kj{;t,^) +K+(-,t,Aj)) (2.34) 



We observe that the leading term in the r.h.s. of (12.311) can be more conveniently written in 
the form 



-f- e-*^^ = ^m-jp^^Vi^M^e-^f (2.35) 
Aj \/5m5E 

n 2 8E , , n,- , 



6 (R-Rj) 



iftiR) = e 5^ +i ir« ) p. = n a ^ p. = p Q (i _ „ #£) (2.37) 



Then it is clear that the leading term has the form of a free evolution of a wave packet which 
starts at t — Tj from the position aj of j th oscillator, with mean momentum Pj. Notice that 
under our assumptions Pj ~ P > 0. 

In particular (I2.3ip gives a precise meaning to the qualitative statement concerning the approx- 
imate behavior of /^(t) made in section 1. 
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3. Second order approximation 



In this section we fix t > r 2 and consider the second order terms corresponding to both oscil- 
lators in some exited states, i.e. terms of the type (see formula (I1.36P ) 

- [ t dsr njni (t-s)v n a ; Q f^ ± (s) 

Jo 

= 4 T n inj (t) [ dse**"' fds'e^'X^is')^ (3.1) 
11 Jo Jo 



i^T ninj (t)h%{t) (3.2) 



for j, I = 1,2, I 7^ j. Proceeding as in lemmas 2.1 and 2.2, a straightforward computation in 
the case a 2 > yields 

hp) = j'ds J'dJ JdH Jdr, G±(-,e,s',?7, S )e^^^ + ' jA; ^^ (3.3) 
Gf t (R, s', V , s) = g 3 (09i(v)e^ (s ' e+sr,2+2s ^ 4>f(R, & s', V , s) (3.4) 
iP±(R, t s>, ry, s) = ^e-^-^T* (3.5) 



R2 = -^-(^' + V s), P£ = Pb±*{£ + V ) (3.6) 

M7 7 

where gj and have been defined in (I2.6I) and (I2.14I) respectively. In the case a 2 < the same 
representation formula 03. 3D holds if we replace A 2 , r 2 with — A 2 , — r 2 , where A 2 = |a 2 |7 _1 , 
r 2 = |a 2 |v c j" 1 . In both cases, we shall discuss the asymptotic behaviour of hi(t) for A 1; A 2 — ► 00. 
The integral (13. 3p contains a rapidly oscillating phase and moreover the phase has exactly one 
critical point. Therefore the behaviour strongly depends on whether or not the critical point 
lies in the integration domain. We shall analyze separately the two cases a 2 > and a 2 < 0. 



3.1. The case a 2 > 0. 

We distinguish the four possible cases: i) h^t), ii) h^t), iii) h^it), iv) hf 2 (t). It is eas- 
ily seen that the point (£0, s' , r] , sq) where the phase is stationary is: (g 2 , r 2 , q\, Ti) for i), 
(— Q2, — t~2, — qi, — Ti) for ii) and (— qi, — n, — 92, — r 2 ) for iii). In all three cases the stationary 
point of the phase does not belong to the domain of integration and then the integral rapidly 
decreases to zero for Ai,A 2 — > 00. On the other hand in the case iv) the stationary point is 
(qi, 7*1, g 2 , r 2 ), i.e. it belongs to the domain of integration and therefore there is a leading term 
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of order (AiA 2 ) 1 which we shall compute. In the next proposition we study the cases iv) 
following the same line of the proof of proposition 2.4. 

Proposition 3.1. For a 2 > and t > r 2 we have 

k+l2{t) = ^^e-^-^GU; qi,n, q 2 , r 2 ) + -L ft+ (•, t, A lf A 2 ) (3.7) 

where TZ^ 2 (R, t, Ai, A 2 ) zs a bounded function of A±, A 2; whose estimate will be given during the 
proof. 

Proof. Let us introduce the change of coordinates x = £— qi, y = s'—T\, w = r\— g 2 , 2 = s— r 2 in 
(13.31) and the shorthand notation G(x, ?/, w, z) = e~ lAiqi ~ tA2q ' 2 Gi 2 (R, x+qx,y+Ti,w + q2, z + r 2 ). 
Then 



t— T2 rz+T-2—Ti 



— T2 J —71 



^•12^) = / °^ / dy I dx I dw G(x,y,w, z)e l T i xy+l T 2 WZ (3. 



a, 6— >oo 



t-72 /-2 + T2-T1 /-a /"ft A A 

1 — ra+i — luz 



G(0, 0,0,0) lim / cfc / dy dx dwe ^ xy+l ^ 



t~2 J —T\ J —a J ~b 



pt-T 2 PZ+T2-T! na pO Ai A2 

+ lim / dz/ dy / dx / dw(G(x,y, 0,0) - G(0, 0,0,0)) e l ^ xy+l ^ wz 

a,b->coJ_ T2 J_ Tl J_ a J_ b 

+ lim / dz / dy dx dw (G(x, y, w, z) — G(x, y, 0, 0)) e* T i xy 1 ^i wz 

a,b^ooJ_ T2 J_ n J_ a J_ b 

= (I) + (//) + (III) (3.9) 
The term (I) is the leading term and it can be easily computed 

Aty 2 tit 2 . f f ~ T2 f z+T2 - T1 sia^-bz sin ^ay 



Att 2 T T ft-Ti PZ+T2-T1 

(I) = — J-*G(0, 0,0,0) lim dz dy 
-G(0, 0,0,0) 



T2 71 

7TZ 7T7y 



47T 2 Tir 2i 

AiA 



2 



47rTiT 2 



e -,A l9l - l A 292G + ( ^ gi)Ti)g2)T2) (310) 
A]A 2 

Concerning the term (II) we have 

27TT9 /"*~ T2 Sin ^2-62 /-^a-ri pa . 

(//) = — - lim / cte 3_ / dy / dx(G(x, y, 0, 0) - G(0, 0, 0, 0)) e n y 



A 2 0,6— ►(» ./— -n. 7T.Z 

72 

A 



-n 

9 7rr /" pTl—Ti A 

""'"/dx/ dy(G(x,y,0,0) - G(0, 0,0,0)) e 1 ^ (3.11) 
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The r.h.s. of (13.111) can be estimated using (12.181) . (12.2 II) and the result is 



I I < 
C 2 



^^G(.,0,0,0)|Ui + \\d x d y G(,0,0,0)\\^ + Jdx\\dld 2 y G(x,;0,0)\\L2 
\ d lGi 2 (R,^,T ll q 2 ,T 2 ) \ + Jd£ \d^d s 'G^ 2 (R,^T U q2,T 2 )\ 

1/2 ! 



A 2 A? 

+ jdi ^lds'\dld 2 s/ Gt 2 (R^,s',q 2 ,T 2 )\ 



(3.12) 



where C 2 is a constant depending on n, t 2 . The term (IiT) can be more conveniently written 

as 



(III) 



dw I dz L(w, z)e T z wz 

Z+T2—TI 

dx I dy (G(x, y, w, z) 



L(w, z) = 

where L(0, 0) = 0. Using (12~19|) . fl2T22l we have 



G(x,y,0,0)) 



(3.13) 
(3.14) 



(J/J) = (IV) + (V) 



(3.15) 



where 



?7r?Y 2 

(/V) = ^4^(0, o) 

|(V-)| < (^||^^(-,0)|Ux-K||^^J^(-,0)|U 1 -H||^^^(-,0)|U 1 -H /c/^H^c/^^C-^, 



(3.16) 



(3.17) 



and C 3 depends on t, r 1; r 2 . Taking into account (13. 14j) we also obtain 



\(v)\<%i\\diG(-,-,-m^HdidM-,-^m^HdidiG(-,,-m^ 



A| 



»t-T 2 



+ dw dx dy\\d w d z G(x,y,w,-)\\L2 



ds' dndr]\d 3 v Gf 2 (R,^s',r],T 2 )\+ / ds' dndri\d 3 v d s Gf 2 (R,^s',r],T 2 



9i 



+ / ds' dUdr]\d v d 2 s Gf 2 (R,^,s',r],r 2 )\+ ds' dudrji ds\d 3 d 3 s Gf 2 (R,£, s' ,r], s) 
Jo J J Jo J J Vo 



1/2 



(3.18) 



Concerning d w d z L(0,0), a straightforward computation gives 
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T2—T1 



-xy 



d w d z L(0, 0) = J dx d w G(x,T 2 — T\, 0,0) e* T i ^ T1 ^ x +JdxJ dy d w d z G{x,y, 0,0) e T i ' 

= (/^) + (/y 2 ) (3.19) 

In (13.191) we integrate by parts in the first integral and use (12. 17ft . (I2.20p in the second integral. 
Then 

< ■^^V||4^G'(-,r 2 -r 1 ,0,0)|Ui 



Al T2 - Tl 



l n 



d£ |d^G| 2 (£,r 2 ,g 2 ,T 2 ) 



Ai r 2 - ri 

\(IV 2 )\ < \\d w d z G(-,0,0,0)\\ L i+ ldx\\d x d y d w d z G(x,;0,0)\\ L * 



(3.20) 



Cx 
Ai 



d£,\d v d s Gf 2 (£,T 1 ,q 2 ,T 2 )\ + J d ^{j ds \d(:d s >d v d s Gf 2 (£,, s' ,q 2 ,r 2 



1/2- 



(3.21) 



and Ci depends on r x , r 2 . Taking into account (I3TT2D . fl3TT8|) . fl3T201 . (l3T2Tj) we get <K7\} . with 
an explicit estimate of 7Zf 2 (R, t, A l5 A 2 ), and this concludes the proof. 



□ 



Let us consider the cases i),ii),iii) where the stationary point of the phase lies out of the 
integration region. In such cases, exploiting repeated integration by parts, one can show that 
(13.31) is O(A^), for any integer k, for Ai — > oo. Since the error term in (13. 7ft is 0(A]~ 3 ), in the 
next proposition we shall limit to k = 3. 

Proposition 3.2. For a 2 > 0, t > r 2 we have 

h a jl {t) = ^n%{;t,K x ,K 2 ) (3.22) 
for a — ±, j — 2,1 — 1 and a = — , j = 1, 1 = 2, where 

\n%(R,tM,k 2 )\<C jdsjds' jdtjdn ( \dlGl{R,t,s\r],s)+d\G%{R,Z,s\ri,s)\^ (3.23) 
and C depends on r 1 ,r 2 . 



Proof. Let us define r = Tl 1 2 T2 and write 
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h+(t) = J\Jds'Jd£ Jd V G+ 1 (-,C,s',v,s)e lA ^^ s ' )+lAld ^ ri ' s) 

+ fds fds'fdi fdr ] GU--X,s' ) r 1) s)e lA ^^ s ' )+lK ^^ s) 
Jtq Jo J J 

= (A) + (B) 

In (A) we integrate by parts three times with respect to the variable £ and we obtain 



(3.24) 



(3.25) 



Integrating by parts with respect to rj in (B) we have 



IB) 



z 3 A? 



< — 



A? (r - nf 



t rs 



ds' d^ \dr){dlGi x {;t,s',r), s))e 



iK 2 6l (e,s')+iAi6>+ (tj,s) 



I ds / ds d^ jdrj Id^Gzii'i £> s \ Vi s ) 

TQ JO J J 



(3.26) 



From (13.25j) and (13. 26ft we get the estimate for h^iit). 
For the estimate of h^t) it is sufficient to notice that 



-«Ai 



_ e iA 2 e-(e,s')^3 e iAi0J-(i7,«) 



(3.27) 



and to integrate by parts three times. The estimate of h 12 (t) is analogous and then the proof 
is complete. 



□ 



Taking into account ( 13.21) . ( 13.71) . ( 13.221) . we obtain for ni,n 2 7^ 0, a 2 > and t > r 2 
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& 2 (t) = + JL K®(.,t, A u A 2 ) (3.28) 

^(2) = _4^2^Z^ e -i(n 1+ n 2+ lM^^ ( 3i29 ) 
AT 

tpi2 = dtki 9i> r i)92, r 2 ) (3.30) 

^(•,t,A 1 ,A 2 ) = z^r nin2 (t) ^ ^(.,t,A l9 A 2 ) (3.31) 

o=± 
j,l=l,2,j^l 

Notice that the leading term in (13.281) can also be written as 



TT^ tK ° ^ = -^ 2 J^Je ^ II ^)(<Mo)fe) e-^t 2 (3.32) 

,. 5E n%5E 5E , . ri\ n 2 

Vviit) = -±— + + n x n 2 — - m + n 2 + l wi + + (3.33) 

2 OTi 2 OT 2 OT 2 OTi OT 2 

Kf (R-R 12 ) 2 , • Pi 2 D 

if)+{R) = -= e R , R 12 = (n iai + n 2 a 2 ) 5E, P 12 = P [l - (m + n 2 ) SE] 

(3.34) 



3.2. The case a 2 < 0. 



Here the two oscillators are on the opposite sides with respect to the origin and one can easily 
check that the point (£ , s' , i] , s ) where the phase in (13. 3p is stationary is: (gi, n, g 2 , — r 2 ) 
for (-gi, -ri, -q 2 ,T 2 ) for (g 2 , -r 2 , g x , r x ) for ^(t), (-g 2 , r 2 , -g l5 -n) for /%(*). 

Since none of these points belongs to the domain of integration we can show that hjj (t) is always 
rapidly decreasing to zero for A l5 A 2 — > oo. 

Proposition 3.3. For a 2 < 0, t > t 2 and any integer k > 2 we have 



hp) = j E Qf l (-,t,A 1 ,A 2 ), j, I = 1,2 (3.35) 



where 



Qfi(R,t,A 1 ,A 2 )\ < JdsJds'JdtJd-n ({d'Gf^ts'iV,^ + {d^iR^s'^s) 



(3.36) 
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Proof. The proof is an immediate consequence of k integration by parts and a trivial estimate. 

□ 



From the above proposition we conclude that for n l5 n 2 ^ 0, a 2 < 0, t > r 2 and any integer 
k > 2 we have 



/£U*) = ^fi C,) 0,*,A 1 ,A 2 ) 
C«(.,t,A 1 ,A 2 )=4r B1B9 (t) J2 G?iO.*,Ai,A a ) 



(3.37) 
(3.38) 



a=± 



4. Joint excitation probability 

We are now in position to compute the joint excitation probability of the two oscillators in 
the two cases a 2 < and a 2 > 0. As a preliminary step, we need a pointwise estimate of the 
derivatives of G^j with respect to the variables £, rj. 
It is convenient to introduce the following notation 





ht 




a 




b 




M7 2 ' 




s 


= to, 


s' 


X 


= <T~ l R, 


z 



tit ° /A 1\ 

c = - (4.1) 



M70" 7 
= t/3, (4.2) 
x+b(f3£ + ari) (4.3) 

We notice that, for t of the same order of magnitude of r 2 , the constants in (14. ip are of order 
one; moreover the rescaled variables a,j3 satisfy < a, (3 < 1. 



Lemma 4.1. We have 



d$Gji(R, C, r/ ; to) I + |dj?G*(i2, £, *ft V, ta) \ 



< c 



k 



* m=0 m=0 



where 



2t fe (t) = (1 + a 2 + b 4 ) fc/2 (1 + b 2 + c 2 ) fe/2 (1 + a) 



(4.4) 



(4.5) 
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Proof. Exploiting the above notation we can write 

G%^^ ± g j g l et' (4.6) 
where ^ = ^(R), 9j = 9i = 9i(v) and 

4> = + \rf + ^Vj-b 2 (^-e + y?? 2 + cW) -bztft + a V )+icx{t + v ) (4.7) 

Let us compute the derivative of order k with respect to rj. 



m—r\ \ / 



m=0 

k 



E (t) ^ E 

m=0 ^ ' n iP 



Til I 



n! p! 2p 

n+2p=m 



A straightforward computation yields 



< [v^Tc*!*! + (a+ bc)(|£| + M) 
= E («) (*' + c 2 )^|zr-«(2a)'(|ei + M) 9 

g=0 ^ ' 

< c(l + b 2 + c 2 ) fe/2 (l + a) k (z) k (O k (v) k (4-9) 
|d 2 0| p < (1 + a 2 + b 4 ) fc/2 (4.10) 



.2 



|^V| = ^:e-T (4.11) 
vo- 



Using gU), (|QD]1 . flCTl in (gSJ we obtain the estimate 



|^G±| < cSU(t)^ <^e"^ (0 fc |ft-(0l fa>*£ K^)l ( 412 ) 

* m=0 



Following exactly the same line we also find the corresponding estimate of and this 

concludes the proof of the lemma. 

□ 



Finally we can prove our main result. 
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Proof of theorem 1. We start with a detailed estimate of V n n (t). Taking into account 
(13371) . (EI38D, fl3T36|) we have 



P» 1W (*)<^E JdR\Q^{R,t, At,A a )\- 

4A 4 /" 
^"^Af E JdR\Q c ; l (R,t,A 1 ,A 2 ) 



a=± 



o </ o 



n 2 



ds ds' d£ d V {\d k v G%\ + \dp%\) 



- 1 StSh'S^ I dv S dR + 



1/2* 



(4.13) 



where in the last line we have interchanged the order of integration and used the Schwartz 
inequality. Exploiting the estimate (I4.4p we find 



A J 



Af " 4 h A Aj 



^■(nA' 2 |y,||^, ; | 



|2 



Jdsjds'{^ JdR (z) 2k e- 2 



1/2 



-i 2 



— A 21, — I 



c A 4 t 



4+4 



Af - 4 £ 4 A 4 



c 



A 



Af " 4 \VSrn6Ej W \fli 



(4.14) 



It remains to evaluate the two norms in (14. 14ft . Recalling the definition of gj(£) (see (I2.6P ) we 
have 



1 9 3 II wb 1 - 



k m 

^EE 

m=0 p=0 

<c||y||rj 



EE 

m=0 p=0 



III 



V J2tv 



l >n j ( Po\\w k ' 1 



d£(0 k C~ p n£K(<Mo)(0 



dx\x\ m - p \V(x)\ <(<Mo)(0 



Inserting (14. 15j) in (I4.14p we finally get the estimate (11.17P with 



K y/Sm5Ej VW W 
Let us consider V+ n2 (t). From fl3T28|) . fl33TD . fl3T32|) . (ET33I) we have 



(4.15) 



,1 (4-16) 
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KJt) = 167r4 ^(TOf) ^ 2 \9M92(Q2)\ 2 + S nin2 (t) (4.17) 
where S nin2 (t) is a correction term of order Aj" 1 . In fact 



\S nin2 (t)\ <j^-^L_M\ gi { qi )g 2 {q 2 )\ (JdR\TZ^(R,t,A 1 ,A 2 )\^j 
+±-JdR\n {2 \R,t, A 1: A 2 )\ 2 

' A " X \^\9M)92{q2)\ (sup [dR\K%(R,t, A U A 2 ^ 

1 ./ 



1/2 



Ai <5m<5£ ft 2 A? 
+^^Z74 SU P I dR\K%(R,t, A U A 2 



Aur*(*) (4.18) 



i /V 4 



Ai \e 



where 



2 / \2 



A/" ^1(51)^(92)1 sup jdR \K%{R, t, A l5 A 2 ; 



vV2 

|2 



+ icO(S?s/' wlwtiA, ' A,) 



(4.19) 



The proof of (11.181) . (11.201) is complete if we notice that the quantity 



^sup [dRlll^t, A U A 2 )\ 2 (4.20) 

t a,j,l J 

can be estimated following the line of the previous case. The explicit computation is straight- 
forward but rather long and tedious and we omit the details. 

□ 

5. Appendix 

Here we give a proof of lemma 2.3 (see e.g. [BHj for analogous computations). 
Proof of lemma 2.3. Let us decompose J(A) in the following form 
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.7(A) = jdx J dyf(x,0)e tAxy + jdx J dy {f{x,y) - f(x,0))e lAxy 

= 4 [dx f(x, 0) e^-e-^ + i 7 ^ y)- 0) ^ 

Aj re AJ J_ u y 

= i(/C n (A) + /C 12 (A)) (5.1) 

where an explicit integration in the first integral and an integration by parts in the second 
integral has been performed. Thus we have (12.171) with /Ci(A) = ^^ =1 /Ciy (A) . The estimate 
of /Ci(A) is easily obtained if we write 

d x f(x,y)-d x f(x,0)=y d6 d x d y f(x,y9) (5.2) 

Jo 

and then use the Schwartz inequality. 

In order to prove (12.181) we reconsider /Cn(A) and /Ci 2 (A). In particular we have 

giA/ix g— iAvx r ^iA/ix ^—iAvx 



£ u (A) = -i/(0,0) /dx z /dx(/(x,0)-/(0,0)) 



XI X 



2,/(o, o) - 1 /('■") -/("■ (£!!! + e - Ara 



A J X 2 \ pi V 

= 27 r/(0,0) + ix: 21 (A) (5.3) 

where we have explicitely computed the first integral and we have integrated by parts in the 
second integral. Concerning /Ci2(A), we observe that it is of the same form as J7"(A) and then we 
can repeat the procedure. Denoting rj(x,y) = dx ^ x ^~ dx ^ x ^ e ^ x v ^ with r](x, 0) = d x d y f(x, 0), 
we obtain 

A J x A J J_ u y 



Ay x Ay ./_„ 

= i(/C 22 (A) + /C 23 (A)) (5.4) 

and the asymptotic formula (12.181) follows, with /C 2 (A) = ^j=i^2j(A). 
The estimate of /C 2 i(A) is obtained if we write 

f{x,0)-f{0,0)-d x f{x,0)x = -x 2 [ d69dlf(x6,0), (5.5) 

Jo 
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the estimate of is trivial and for /C23(A) it is sufficient to notice that 

d 2 J(x,y) - d 2 J(x,0) - d 2 x d y f(x,0)y = y 2 I ' d66 [ d( d 2 x d 2 y f(x,y90 

Jo Jo 

and to use the Schwartz inequality. Then the estimate (I2.2ip for /C 2 (A) is proved. 
Finally we shall prove (12.191) . An integration by parts in /C2i(A) yields 



(5.6) 



r ( w 21 L fM-f(0,0)-d x f(x,0)x + d 2 J(x,0)^fe iA ^ e~^\ 
/C 21 (A) = - jdx ^— —) 

= i/C 3 i(A) (5.7) 
For /022(A) we proceed as in (15.31) and we obtain 

giA/ix g— iAvx p ^iAfix ^—iAvx 



dx + jdx (d x d y f(x,0)-d x d y f '(0,0))' 

. . i f 7 d x d v f(x,0) -d x d v f (0,0) -d 2 d v f(x,0)x f e iA " x e~ iAux 
= 2irid x d v f (0,0) - - dx yjy ' ; ^ v ' ; x yJ v - ; + 



A J x 2 \ fj, v 

= 2md x d y f (0,0) + ^/C 3 2 (A) (5.8) 
The last term /C 2 3(A) has the same form as J7"(A) and then following the same argument we get 

„■ [■ _iA.ua; c —iAvx 



/C 23 (A) = ^ Jdxd 2 x d 2 y f(x,0)- 



x 

2 



_i_ f dx r dy 4f(x,y) - dlfjxfi) - d 3 x d y f(x,0)y - d^fjxfi)^ ^ 
A J J- v y 3 

= i(/C33(A)+/C 3 4(A)) (5.9) 
and fl2TT9|) is proved with /C 3 (A) = Ej=i ^3j( A )- 

The estimate (I2.22p for ^(A) is easily obtained following the same line as before and we omit 
the details. 

□ 
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